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ABSTRACT 

We report the detection, from observations using the James Clerk Maxwell Telescope, of CO 
J = 3 — > 2 transition lines in six carbon stars, selected as members of the Galactic Halo and 
having similar infrared colors. Just one Halo star had been detected in CO before this work. 
Infrared observations show that these stars are red (J-K >3), due to the presence of large dusty 
circumstellar envelopes. Radiative transfer models indicates that these stars are losing mass 
with rather large dust mass-loss rates in the range 1-3.3 x 10 _8 M Q yr _1 , similar to what can 
be observed in the Galactic disc. We show that two of these stars are effectively in the Halo, 
one is likely linked to the stream of the Sagittarius Dwarf Spheroidal galaxy (Sgr dSph), and 
the other three stars certainly belong to the thick disc. The wind expansion velocities of the 
observed stars are low compared to carbon stars in the thin disc and are lower for the stars in 
the Halo and the Sgr dSph stream than in the thick disc. We discuss the possibility that the 
low expansion velocities result from the low metallicity of the Halo carbon stars. This implies 
that metal-poor carbon stars lose mass at a rate similar to metal-rich carbon stars, but with 
lower expansion velocities, as predicted by recent theoretical models. This result implies that 
the current estimates of mass-loss rates from carbon stars in Local Group galaxies will have 
to be reconsidered. 

Key words: circumstellar matter - infrared: stars — carbon stars — AGB stars — stars: mass 
loss 



1 INTRODUCTION 

Stars with initial masses in the range 0.8-8 Mq end their life with 
a phase of catastrophic mass-loss. During the asymptotic giant 
branch (AGB) phase, they develop a superwind leading to mass- 
loss rates up to 1O _4 M0 yr _1 . This superwind enriches the ISM 
with newly synthesized elements. 

The mass-loss mechanism of AGB stars is likely due to pul- 
sations from the star and radiation pressure on dust grains. Shocks 
due to pulsation extend the atmosphere, so that the material ejected 
by the star becomes dense and cold enough for dust to form. Due 
to its opacity, dust absorbs the radiation from the star and is driven 



* Based on observations made at the James Clerk Maxwell Telescope un- 
der program M08AU01 

f Based on observations made at the Very Large Telescope at Paranal Ob- 
servatory under the program 082.D-0836 
| E-mail: eric.lagadec@manchester.ac.uk 



away by radiation pressure, carrying the gas along through fric- 
tion. Theoretical models (Winters et al. 2000), show that the mass 
loss evolves from a pulsation driven regime characterized by a 
low mass-loss rate and a slow expansion velocity to a dust-driven 
regime with a high mass-loss rate and a high expansion velocity 
OSkm.s" 1 ). 

Studying the effect of metallicity on the mass-loss process is 
important to understand the formation of dust around AGB stars 
in the early Universe. At low metallicity, less seeds are present for 
dust formation, so one might expect dust formation to be less effi- 
cient and thus the mass-loss rates to be lower. 

Theoretical work by Bowen & Willson (1991) predicts that 
for metallicities below [Fe/H]= —1 dust-driven winds fail, and the 
wind must stay pulsation-driven. However, observational evidence 
for any metallicity dependence is still very limited (Zijlstra 2004). 
More recent observational (Groenewegen et al. 2007) and theoreti- 
cal works (Wachter et al. 2008; Mattsson et al. 2008) indicate that 
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the mass-loss rates from carbon stars in metal-poor environments 
are similar to our Galaxy. 

To obtain the first observational evidence on mass-loss rates 
at low metallicity, we have carried out several surveys with the 
Spitzer Space Telescope of stars in nearby dwarf galaxies. These 
show significant mass-loss rates down to Z=l/25 Zq (Lagadec et 
al. 2007b; Matsuura et al. 2007; Sloan et al. 2009), but only for 
carbon-rich stars. The current evidence indicates that oxygen-rich 
stars have lower mass-loss rates at lower metallicities. For carbon 
stars, no evidence for a dependency of mass-loss rate on metallicity 
has yet been uncovered. Consequently, (Lagadec & Zijlstra 2008) 
have proposed that the carbon-rich dust plays an important role in 
triggering the AGB superwind. 

The main uncertainty arises from the unknown expansion ve- 
locity. This parameter is needed to convert the density distribution 
to a mass-loss rate. The expansion velocity is also in itself a pow- 
erful tool. Hydrodynamical simulations (Winters et al. 2000), have 
shown that radiation-driven winds have expansion velocity in ex- 
cess of 5kms -1 , while pulsation-driven winds are slower. 

There is some evidence that expansion velocities are lower at 
low metallicity, from measurement of OH masers (Marshall et al. 
2004). However, OH masers are found only in oxygen-rich stars, 
which appear to have suppressed mass-loss at low metallicities. 
We lack equivalent measurements for metal-poor carbon-rich stars, 
which do reach substantial mass-loss rates. For these, the only 
available velocity tracer is CO. Currently, extra-galactic stars are 
too distant for CO measurements. However, there are a number of 
carbon stars in the Galactic Halo, which are believed to have sim- 
ilarly low metallicity. One CO measurement exists for a Galactic 
Halo star: its expansion velocity has been estimated to be ~3.2 
kms -1 through the 12 CO J = 2 — > 1 transition (Groenewegen 
et al. 1997). This is much lower than typical expansion velocities, 
which are in the range 10-40 km s _1 for stars with similar colours, 
and thus optical depths (FigQ}. 

A large number of metal-poor carbon stars have recently been 
discovered in the Galactic Halo (Totten & Irwin 1998; Mauron et 
al. 2004, 2005, 2007). These stars may have been stripped from the 
Sagittarius Dwarf Spheroidal galaxy (Sgr dSph), which has a metal- 
licity of [Fe/H]~ —1.1 (Van de Bergh 2000). These stars are the 
closest metal-poor carbon stars known, and they are bright enough 
to be detected in CO using ground-based millimeter telescopes. 

We have carried out observations of six Halo carbon stars in 
the CO J = 3 — > 2 transition. Here, we report the results of these 
observations. 



2 SAMPLE SELECTION 

Mauron et al. (2004, 2005, 2007) and Mauron (2008) have discov- 
ered ~ 100 carbon stars in the Galactic Halo, adding to the sample 
of ~ 50 described by Totten & Irwin (1998). All of these stars are 
spectroscopically confirmed carbon stars. Only the brightest can 
be detected in the sub-millimeter range. We selected the six stars 
with the highest IRAS 12/im flux observable with the James Clerk 
Maxwell Telescope (JCMT, Mauna Kea, Hawaii). The emission 
from an AGB star at 12^im is due to thermal emission from the 
dust in the envelope. Thus one expects the stars with the largest 
12^m flux to be the brightest in CO. All the observed stars have 
3< J ~ K <4 and have thus a similar dust optical depth and a 
large circumstellar dusty envelope. 

TableQ] describes the observational properties of the selected 
stars. One of these stars, IRAS 12560+1656 (Beichman et al. 1990), 
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Figure 1. Distribution of the expansion velocity for the observed Halo car- 
bon stars compared with stars from the disc with similar J— K colors. 



is the only Halo carbon star that had previously been detected in CO 
J= 2 — > 1 (by Groenewegen et al. 1997). 



3 OBSERVATIONS AND DATA REDUCTION 

We observed the CO J = 3^ 2 (345 GHz) line emission from the 
six selected carbon stars from the JCMT using the heterodyne 
focal-plane array receiver (HARP; Smith et al. 2008). Observations 
were carried with a bandwidth of 1 GHz and a frequency resolution 
of 0.977 MHz, giving a velocity coverage of ±400 kms^ 1 . To in- 
crease the signal-to-noise ratio of our observations, we used a 30" 
beam switch to keep one detector on the source at all times. We 
used the Starlink software to reduce the resulting data. Fig.[2]shows 
the six CO 3-2 lines we detected. 

We obtained mid-infrared photometric observations of two 
stars from our sample using VISIR (Lagage et al. 2004) on the VLT. 
The observations were taken using four filters centred at 8.59/im 
(PAH1, AA=0.42^m), 11.25^m (PAH2, AA=0.59^m), 13.04/an 
(Nell 2, AA=0.22^m) and 17.65/im (Ql, AA=0.83^m). We used 
the standard mid-infrared chop-and-nod technique to remove the 
background emission from the telescope and sky. Standard stars 
were observed just after each observation for flux calibration. The 
data were reduced using the IDL routines developed and described 
by Lagadec et al. (2008). 



4 CO LINE PROPERTIES 

Fig.|2]shows that the six lines we detected have parabolic profiles, 
which arise when the envelopes of the observed stars are spheri- 
cal, unresolved by the telescope, have a constant expansion veloc- 
ity V^exp, and the CO 3-2 line is optically thick (Knapp & Morris 
1985). This can be described as: 



TX(V) = TX (peak) 



1 



V- Vc 
Vn 



(1) 
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Figure 2. JCMT CO 3-2 line of the 6 observed stars. The dashed line represent the best model fit to the data. The velocities in abscissa are Vlsr 
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Table 1. Observed Halo stars targets: names, adopted coordinates, photometry and distance. J, H, K s are taken from 2MASS. fg.59, fi3.04 an d fi7.65 are tne 
VISIR/VLT flux at 8.59, 13.04 and 17.65 fim respectively. fi2 is the IRAS 12|tm flux. Z is the distance above the Galactic plane. 



IRAS name 


Other name 


RA 


Dec 




K 


J-K 


f8.59 


fl3.04 


fl7.65 


fl2 


[12] 


Z 






(J2000) 




mag 


mag 


Jy 


Jy 


Jy 


Jy 


mag 


kpc 


IRAS 04188+0122 


COS 6075 


04 21 27.25 


+01 


29 13.4 


6.420 


3.280 


4.99 


2.67 


1.64 


3.37 


2.31 


-2.63 


IRAS 08427+0338 


CGS 6306 


08 45 22.27 


+03 


27 11.2 


6.255 


3.410 








6.50 


1.60 


2.03 


IRAS 11308-1020 


CGS 3052 


1 1 33 24.57 


-10 


36 58.6 


4.568 


3.830 








57.37 


-0.77 


1.41 


IRAS 16339-0317 


CGS 3716 


16 36 31.70 


-03 


23 37.5 


6.098 


3.906 


12.49 


10.83 


5.60 


14.57 


0.72 


1.73 


IRAS 12560+1656 


CGS 6500 


12 58 33.50 


+ 16 


40 12.0 


7.820 


3.480 








0.77 


3.91 


8.96 


IRAS 18120+4530 




18 13 29.6 


+45 


31 17.0 


6.710 


3.809 








7.86 


1.39 


2.20 



Table 2. CO J=3-2 lines properties of the observed Halo stars. Typical errors on the expansion velocities estimates are 0.1 km.s 1 for the stars with expansion 
velocities above 10km. s~ 1 and 0.2s _1 for the other ones 



Adopted name 


peak 




^cxp 


FWHM 


flux 


rms 


Integration time 




K 


km.s 


km.s 


km.s - 1 


K. km.s- 1 


K 


mn 


IRAS 04188+0122 


0.04 


4.8 


11.5 


17.3 


0.77 


0.009 


616 


IRAS 08427+0338 


0.03 


20.9 


16.5 


17.4 


0.55 


0.010 


45 


IRAS 11308-1020 


0.79 


20.8 


11.5 


12.9 


10.86 


0.049 


90 


IRAS 16339-0317 


0.19 


-83.1 


8.5 


12.6 


2.59 


0.030 


182 


IRAS 12560+1656 


0.03 


88.6 


3. 


4.5 


0.14 


0.008 


453 


IRAS 18120+4530 


0.08 


-268.0 


6.5 


8.3 


0.72 


0.019 


227 



where T\ (peak) is the peak emission temperature, Vo is the expan- 
sion velocity of the envelope, and Vc the velocity of the line. We 
fitted this function to the observed data with a least-squares tech- 
nique. Table[2] lists the resulting properties of the observed lines. 
The CO lines detected for IRAS 12560+1656 is rather weak, but 
can be considered as a real detection as it was already detected in 
another transition by Groenewegen et al. (1997). 

FigO shows the distribution of expansion velocities for our 
sample. As expected, the flux measured in the CO 3-2 line corre- 
lates strongly with the IRAS 12/im flux (Fig.[4}. The best linear fit 
to the data is obtained for: 

Fco = 0.194 ± 0.009 x f 12 - 0.314 ± 0.214 (2) 



5 DUST MASS-LOSS RATES 

Two techniques are usually used to estimate dust mass-loss rates 
from AGB stars: by measuring the near- and mid-infrared flux of 
the stars (Whitelock et al. 1994, 2006; Lagadec et al. 2008) or by 
using radiative transfer models. We apply both methods here. We 
estimated the mass-loss rates by measuring the K s — [12] colour, 
where K s is the 2MASS 2.16^im magnitude and [12] the IRAS 
12/^m magnitude, according to the relation from Whitelock et al. 
(2006): 

log(Mtotai) = -7.668 + 0.7305(A" - [12]) 
-5.398 x W- 2 {K - [12]) 2 
+1.343 x W~ 3 {K - [12]) 3 , (3) 

where [12] is the IRAS 12/^m magnitude assuming a zero- 
magnitude flux of 28.3 Jy. 

We obtained an alternative estimate of the dust mass-loss rates 



using a radiative transfer model. We used the radiative transfer 
code DUSTY (Ivezic & Elitzur 1997). This code solves the 1- 
dimensional problem of radiation transport in a dusty environment. 
For all of our models, we assume that the irradiation comes from a 
point source (the central star) at the centre of a spherical dusty en- 
velope. The circumstellar envelope is filled with material from a ra- 
diatively driven wind. All of the stars are carbon-rich, and the dust 
consists of amorphous carbon. We did not add SiC to our model, 
except for IRAS 04188+0122 and IRAS 16339-0317, for which a 
fit was impossible without SiC. Optical properties for these dust 
grains are taken from Hanner (1988) and Pegourie (1988) for amor- 
phous carbon and SiC, respectively. The grain size distribution is 
taken as a typical MRN distribution, with a grain size a varying 
from 0.0005 to 0.25/im distributed according to a power law with 
n(a)<xa~ q with q=3.5 (Mathis et al. 1977). The outer radius of the 
dust shell was set to 10 3 times the inner radius; this parameter has 
a negligible effect on our models. 

To model the emission from the central star, we used a hydro- 
static model including molecular opacities (Loidl et al. 2001 ; Groe- 
newegen et al. 2007). Our aim was to fit the spectral energy distri- 
bution, defined by the 2MASS photometry (J, H and K), our VISIR 
mid-infrared data, and the IRAS photometry (at 12 and 25/im) to 
estimate the dust mass-loss rates. We fixed the dust temperature at 
inner radius (1220K), the mass ratio of SiC to amorphous (0 %) 
carbon dust and the effective temperature of the central star, T e ff 
(2800K), unless a satisfactory fit could not be obtained with these 
parameters. 

DUSTY gives total (gas+dust) mass-loss rates assuming a gas- 
to-dust ratio of 200. Estimating the gas-to-dust ratio of these stars 
is beyond the scope of this paper and will be the subject of a forth- 
coming paper using CO observations of these stars in other transi- 
tions and Spitzer spectra. We thus prefer to estimate dust mass-loss 
rates. The dust mass-loss rate can be obtained by dividing by the 
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Table 3. Dust mass-loss rates for the observed stars and parameters obtained from our DUSTY models. M co \ and Af dusty are the dust mass-loss rates (in 
MQyr" -1 ) from Lagadec et al. (2008) and our DUSTY models, respectively. The distance we adopted is the average between the distances Di and D2 described 
in Sectionl6!2l 



Target Luminosity MduSTY ^-eff ^in SiC/AMC r (0.55/^m) M co \ Distance 

(L ) (lO-SMayr- 1 ) (K) (K) (lO^Meyr- 1 ) kpc 
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2800 


1200 


0.1 


10. 


1.6 


6.2 
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2800 
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10. 
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Figure 4. The linear correlation between the IRAS 12/im flux and the ob- 
served CO J=3-2 flux. 



Figure 5. Total (gas+dust) mass-loss rate as a function of the absolute bolo- 
metric magnitude. 



gas-to-dust mass ratio, assuming that the gas and dust expansion 
velocities are the same. Here we prefer to talk about dust mass-loss 
rates as inputs from our models as we are fitting infrared colours, 
and thus dust emission. The expansion velocity is an output for the 
DUSTY models, and we scaled our results using the expansion ve- 
locities we measured. Table[3]gives the results of the fits. The IRAS 
12^tm flux is an outlier for IRAS 16339—0317. This is certainly due 
to variability or to the fact that the aperture of the IRAS satellite is 
larger than the VISIR one. 

6 DISCUSSION 

6.1 Mass loss and evolution 

The detection of resolved CO lines in our six targets indicate that 
all of these stars are losing mass. Our VISIR observations, as well 
as IRAS and 2MASS observations, show that these stars are very 
red, which indicates the presence of circumstellar dust. To study 
the impact of the mass loss on the evolution of the studied stars, we 
compared the measured mass-loss rates with nuclear reaction rates 



and the classical single-scattering limit (see e.g., van Loon et al. 
1999 and Lagadec et al. 2008). 

Fig[5] shows the total (gas+dust) mass-loss rate as a function 
of the absolute magnitude. The absolute magnitude is taken from 
our DUSTY models. The total mass-loss rate is also an output from 
our DUSTY models, assuming a typical gas-to-dust mass ratio of 
200. The dotted line shows the rate at which mass is consumed by 
nuclear burning and the dashed line the classical single-scattering 
limit as described by Lagadec et al. (2008). This shows that no 
multiple scattering is needed to explain the observed dust mass- 
loss rates. The mass-loss rates we measure are also well above 
the nuclear burning rate line, indicating that the evolution of these 
stars will be governed by the mass-loss process. Similar conclu- 
sions were drawn from observations of metal-poor AGB stars in 
the LMC, Fornax and Sgr dSph (van Loon et al. 1999, Lagadec et 
al. 2008). 

6.2 Distances 

Most of the methods to determine mass-loss rates rely on an ac- 
curate distance to the observed stars. To determine these distances, 
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Figure 3. Radiative transfer DUSTY models of the spectral energy distributions of the observed stars. The dashed lines represent our best models, and the 
symbols are flux from 2MASS, IRAS (12 and 25/im) and our VISIR observations. Filled squares, open squares and triangles represent 2MASS, IRAS and 
VISIR/VLT data, respectively. 
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Table 4. Galactic location of our observed stars. Di, D2 and D3 are the distance to the sun using methods described by Sloan et al. (2008), Mauron et 
al. (2008) and the period-luminosity relation respectively. This last method was not applied for three stars without known periods from the Northern Sky 
Variability Survey (NSVS, Wozniak et al. 2004). 



Adopted name 


1 


b 


Period 


Di 


D 2 


D3 


<D> 








days 


kpc 


kpc 


kpc 


kpc 


IRAS 04188+0122 


192.1775 


-31.9867 


359 


7.3 


6.0 


6.4 


6.5±0.6 


IRAS 08427+0338 


223.4859 


+26.8173 


288 


6.6 


5.3 


5.1 


5.5±0.8 


IRAS 11308-1020 


273.6969 


+47.7772 




2.8 


2.1 




2.5±0.5 


IRAS 16339-0317 


012.7346 


+27.7944 




5.6 


4.2 




4.9±1.0 


IRAS 12560+1656 


312.2528 


+79.4127 




13.3 


10.7 




12.0±1.9 


IRAS 18120+4530 


073.0530 


+25.3482 


408 


7.6 
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6.7 
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fitted line, compared to 0.5 1 magnitudes for the SMC sample. The 
two fitted lines yield nearly identical distances at J— K = 2, but as 
the colour grows redder, the samples diverge from each other. At 
J— K = 5, the difference in Mk is a full magnitude. For the colours 
in our sample, the two methods yield results differening by 0.41 to 
0.64 magnitudes, comparable to the spread in the SMC sample. The 
different slopes in the two samples may result from their different 
metallicities, but that is only speculation on our part. 

The third method utilises the period-luminosity relation for 
carbon Miras described by Feast et al. (2006): 



-2.54 log P + 2.06, 



(6) 



With an uncertainty of 0.24 magnitudes. We derived the bolometric 
magnitudes using the equation for bolometric correction derived by 
Whitelock et al. (2006), after converting all of the photometry to the 
SAAO system as described by Lagadec et al. (2008): 



Figure 6. Absolute K magnitude (Mk ) of carbon stars as a fonction of 
their J-K infrared colours. Diamonds and squares represent LMC and SMC 
carbon stars respectively. The vertical dashed line represent the limit under 
which relation Eq|5]is no longer valid. 

we applied three methods. Two are based on near-infrared colours; 
the third is the period-luminosity relationship. 

The first method uses the infrared colours of the observed 
stars, using the relation between Mk and J— K determined by Sloan 
et al. (2008). They found from a sample of carbon stars in the Small 
Magellanic Cloud (SMC) that: 

M K = -9.18 + 0.395(J - #) (4) 

Mauron et al. (2008) compared the near-infrared 2MASS pho- 
tometry of stars in the Halo and the Large Magellanic Cloud (LMC) 
and found a different relationship from the one above between J— K 
and Mk- Their relation gives fainter Mk values at a given J— K 
colour than that of Sloan et al., with the difference increasing with 
redder colours. We have recalibrated the relation for the LMC, us- 
ing the samples of stars confirmed to be carbon-rich with the In- 
frared Spectrograph on Spitzer by Zijlstra et al. (2006), Buchanan 
et al. (2006), Leisenring et al. (2008), and Sloan et al. (2008). We 
find that: 

M K = -9.94 + 0.754(J-A"), (5) 

for J— K colours greater than 2.0. This calibration of the J— K re- 
lation is our second method. Fig. [6] compares the SMC and LMC 
calibrations of the relation between Ma- and J— K. 

The scatter in the LMC sample is 0.81 magnitudes about the 



BC K = + 0.972 + 2.9292 x ( J - K) - 1.1144 x ( J - K) 2 
+0.1595 x (J - Kf - 9.5689 10 3 (J - Kf (7) 

Table|4]presents the obtained distances, Di , D2 and Da respec- 
tively. For those stars without periods, the final estimated distance 
is the average of Di and D2. For those stars with periods, the Di 
and D2 values bracket D3 in two of the three cases. Consequently, 
we first averaged Di and D2, then averaged the result with D3 to 
arrive at our final estimate of the distance. The uncertainties in the 
final estimated distances are the the standard deviation of the indi- 
vidual distances. 

6.3 Galactic location of the stars 

Knowing the distance to the stars and their Vj sr allows us to study 
their location in the Galaxy. All six stars examined here have been 
classified previously as members of the Halo, based on their dis- 
tances from the Galactic plane. Our CO observations give us veloc- 
ity information for these stars, which we can compare to Galactic 
rotation models. Table|4]lists the Galactic coordinates 1 and b. 

Fig[7j shows the location of our six carbon stars on an Aitoff 
projection. The dashed line schematically represents the Sgr dSph 
orbit (Ibata et al. 2001). One of our stars, IRAS 12560+1656, 
lies very close to this orbit. Its VJ sr is in the range of the ob- 
served Vim for stars in the Sgr dSph stream. IRAS 12560+1656 
very likely belongs to this stream. Our observations are thus cer- 
tainly the first CO observations of an extragalactic AGB star. 
Two other stars, IRAS 16339-0317 and IRAS 18120+4530, have 
large negative Vi sr , fully consistent with membership in the Halo. 



8 E. Lagadec et al. 



Finally, IRAS 04188+0122, IRAS 08427+0338, and IRAS 11308- 
1020 have a distance, location and Vj sr consistent with member- 
ship in the thick disc. These three last stars thus have a metallicity 
between that of the thin Galactic disc and the Galactic Halo, the 
average metallicity of the thin disc being ~-0.17 while the one of 
the thick disc is ~-0.48 (Soubiran et al., 2003). Our sample thus 
contains three metal-poor AGB stars and three AGB stars with in- 
termediate metallicity. 

6.4 Low expansion velocities in the Halo 

The present observations allow us to directly measure the expan- 
sion velocity for a sample of carbon stars in the Halo. The stars we 
observed are quite red (3<J— K<4), and have substantial circum- 
stellar dusty envelopes responsible for the observed reddening. Our 
measured expansion velocities are in the range 3-16.5 km s _1 . The 
expansion velocity of carbon stars increases during the evolution 
on the AGB (Schoi'er, 2007), i.e. when the dusty envelope becomes 
optically thicker. To compare the expansion velocities we measured 
in Halo carbon stars with carbon stars in the disc, we took a sam- 
ple of carbon stars in the disc with colours similar to our sample. 
We selected all of the carbon stars with 3<I— K<4 in the exten- 
sive catalogue of CO observations of evolved stars by Loup et al. 
(1993). Fig[T] compares the distribution of expansion velocities in 
the two samples. The Halo carbon stars clearly have a lower mean 
expansion velocity. 

The three stars in the Halo and the Sgr dSph stream have V exp 
in the range 3-8.5 km s , while the three stars associated with the 
thick disc have velocities ranging from 11.5 to 16.5 kms -1 . The 
latter range is at the low end of expansion velocities for AGB stars 
with similar near-infrared colours in the thin disc. 



6.5 Origin of the low expansion velocities 

Section 1631 and [641 have shown that the three stars we observed 
in the Halo and Sgr dSph stream have low expansion velocities. 
The stars we observed in the thick disc have expansion velocities 
intermediate between those typically observed in the Halo and the 
thin disc. This difference could arise from differences in metallicity, 
with the more metal-poor carbon stars having the slower winds. 

Mattsson et al. (2008) and Wachter et al. (2008) have recently 
conducted theoretical investigations of the winds from metal-poor 
carbon stars. Both studies show that metal-poor carbon stars can 
develop high mass loss rates, leading to the formation of a large 
dusty envelope, in agreement with spectroscopic observations from 
Spitzer of AGB stars in metal-poor galaxies (Zijlstra et al. 2006; 
Sloan et al. 2006; Groenewegen et al. 2007; Lagadec et al. 2007; 
Matsuura et al. 2007; Leisenring et al. 2008; Lagadec et al. 2009; 
Sloan et al. 2009). 

Wachter et al. (2008) predict that the outflow velocities from 
carbon stars should be lower in metal-poor environments, because 
of the lower gas-to-dust mass ratio and because the formation of 
less dust leads to less efficient acceleration of the wind outside of 
the sonic region. This interpretation is consistent with our interpre- 
tation that the low expansion velocities we have observed in the 
Halo are due to their low metallicities. 

Hydrodynamical models (Winters et al. 2000, Wachter et al. 
2008) distinguish two types of models. Model A applies to cases 
where the radiation pressure on dust is efficient. Mass-loss rates can 
exceed 10 _7 Mq yr _1 , and expansion velocities can climb above 
~5 kms -1 . Model B applies to cases where pulsations drive the 



mass loss. In these cases, the mass-loss rates and expansion veloc- 
ities are smaller. The mass loss occurs in a two-step process, with 
stars first losing mass due to pulsations, followed by acceleration 
due to radiation pressure on the dust grains. The high mass-loss rate 
and low expansion velocity of TI 32 does not fit either of these mod- 
els, possibly because metal-poor carbon stars can develop strong 
mass-loss from pulsation alone. 



7 CONCLUSIONS AND PERSPECTIVES 

We have detected the CO J = 3— > 2 in six carbon stars selected as 
Halo stars. Only one carbon star had been detected in CO previ- 
ously. Comparison of the infrared observations and radiative trans- 
fer models indicates that these stars are losing mass and produc- 
ing dust. Their mass-loss rates are larger than their nuclear burning 
rates, so their final evolution will be driven by this mass-loss phe- 
nomenon. 

We show that three of the observed stars are certainly mem- 
bers of the thick disc, while one is in the Sgr dSph stream and 
two are in the Halo. The CO observation of the Sgr dSph stream 
star are thus the first identified millimetre observations of an ex- 
tragalactic AGB star. The expansion velocity we determined from 
our CO observations are lower than those of carbon stars in the thin 
disc with similar near-infrared colours. The observed carbon stars 
with the lowest expansion velocities are Halo or Sgr dSph stream 
carbon stars. There is a strong indication that the expansion winds 
are lower in metal-poor environments, which agrees with recent 
theoretical models (Wachter et al. 2008). 

So far, the effect of metallicity on the mass-loss from carbon- 
rich AGB stars has studied primarily from infrared observations 
of AGB stars in Local Group galaxies, mostly with the Spitzer 
Space Telescope. Infrared observations measure the infrared ex- 
cess, which can be converted to a mass-loss rates assuming an 
expansion velocity for the circumstellar material. So far all of the 
mass-loss rates have been estimated using the assumptions that the 
expansion velocity is independent of the metallicity. The results 
presented here show that this assumption needs to be reconsidered. 

We have recently obtained spectra of some of the present 
sample with the Infrared Spectrograph on Spitzer. Combining our 
CO observations and these spectra together and comparing them 
to Spitzer spectra from carbon stars in other galaxies in the Lo- 
cal Group will allow us to quantitatively study the mass-loss from 
carbon-rich AGB stars in the Local Group and its dependence on 
metallicity. This will be the subject of a forthcoming paper. 
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